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£N) , In a two-flavor color superconductor, the SU(3) C gauge symmetry is spontaneously broken by 

diquark condensation. The Nambu-Goldstone excitations of the diquark condensate mix with the 
gluons associated with the broken generators of the original gauge group. It is shown how one can 
decouple these modes with a particular choice of 't Hooft gauge. We then explicitly compute the 
spectral density for transverse and longitudinal gluons of adjoint color 8. The Nambu-Goldstone 
excitations give rise to a singularity in the real part of the longitudinal gluon self-energy. This leads 
to a vanishing gluon spectral density for energies and momenta located on the dispersion branch of 

£SJ ' the Nambu-Goldstone excitations. 

7— I \ 

CN 1 I. INTRODUCTION 

> : 

o : o 

Cold, dense quark matter is a color superconductor ]1J]. For two massless quark flavors (say, up and down), Cooper 

■ pairs with total spin zero condense in the color-antitriplet, flavor-singlet channel. In this so-called two-flavor color 
superconductor, the SU(3) C gauge symmetry is spontaneously broken to SU(2) C ||. If we choose to orient the (anti-) 
color charge of the Cooper pair along the (anti-) blue direction in color space, only red and green quarks form Cooper 
pairs, while blue quarks remain unpaired. Then, the three generators T\, T2, and T3 of the original SU(3) C gauge 
group form the generators of the residual SU(2) C symmetry. The remaining five generators T4, . .. ,Tg are broken. 

& ' (More precisely, the last broken generator is a combination of T§ and the generator 1 of the global U(l) symmetry of 
1 I baryon number conservation, for details see Ref. Q and below). 

According to Goldstone's theorem, this pattern of symmetry breaking gives rise to five massless bosons, the so- 
called Nambu-Goldstone bosons, corresponding to the five broken generators of SU(3) C . Physically, these massless 

■ bosons correspond to fluctuations of the order parameter, in our case the diquark condensate, in directions in color- 
flavor space where the effective potential is flat. For gauge theories (where the local gauge symmetry cannot truly 
be spontaneously broken), these bosons are "eaten" by the gauge bosons corresponding to the broken generators of 
the original gauge group, i.e., in our case the gluons with adjoint colors a — 4, ... ,8. They give rise to a longitudinal 
degree of freedom for these gauge bosons. The appearance of a longitudinal degree of freedom is commonly a sign 

- - 1 that the gauge boson becomes massive. 

In a dense (or hot) medium, however, even without spontaneous breaking of the gauge symmetry the gauge bosons 
already have a longitudinal degree of freedom, the so-called plasmon mode ^J] . Its appearance is related to the presence 
of gapless charged quasiparticles. Both transverse and longitudinal modes exhibit a mass gap, i.e., the gluon energy 
Po — ► nig > for momenta p — ► 0. In quark matter with Nf massless quark flavors at zero temperature T = 0, the 
gluon mass parameter (squared) is 



X 



where g is the QCD coupling constant and /i is the quark chemical potential. 

It is a priori unclear how the Nambu-Goldstone bosons interact with these longitudinal gluon modes. In particular, 
it is of interest to know whether coupling terms between these modes exist and, if yes, whether these terms can be 
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eliminated by a suitable choice of ('t Hooft) gauge. The aim of the present work is to address these questions. We 
shall show that the answer to both questions is "yes" . We shall then demonstrate by focussing on the gluon of adjoint 
color 8, how the Nambu-Goldstone mode affects the spectral density of the longitudinal gluon. 

Our work is partially based on and motivated by previous studies of gluons in a two- flavor color superconductor . 
The gluon self-energy and the resulting spectral properties have been discussed in Ref. [Q. In that paper, however, 
the fluctuations of the diquark condensate have been neglected. Consequently, the longitudinal degrees of freedom of 
the gluons corresponding to the broken generators of SU (3) c have not been treated correctly. The gluon polarization 
tensor was no longer explicitly transverse (a transverse polarization tensor IP" obeys P M IT r " / = IP"' P v = 0), and 
it did not satisfy the Slavnov- Taylor identity. As a consequence, the plasmon mode exhibited a certain peculiar 
behavior in the low- momentum limit, which cannot be physical (cf. Fig. 5 (a) of Ref. It was already realized in 
Ref. |7j that the reason for this unphysical behavior is the fact that the mixing of the gluon with the excitations of the 
condensate was neglected. It was moreover suggested in Ref. j7j that proper inclusion of this mixing would amend the 
shortcomings of the previous analysis. The aim of the present work is to follow this suggestion and thus to correct the 
results of Ref. with respect to the longitudinal gluon. Note that in Ref. (5| fluctuations of the color-superconducting 
condensate were taken into account in the calculation of the gluon polarization tensor. As a consequence, the latter is 
explicitly transverse. However, the analysis was done in the vacuum, at fi = 0, not at (asymptotically) large chemical 
potential. 

The outline of the present work is as follows. In Section || we derive th e transverse and longitudinal gluon 



propagators including fluctuations of the diquark condensate. In Section III we use the resulting expressions to 
compute the spectral density for the gluon of adjoint color 8. Section IV concludes this work with a summary of our 
results. 

Our units are h = c = feg = 1. The metric tensor is g^ v = diag (+,—,—, —). We denote 4-vectors in energy- 
momentum space by capital letters, — (fco,k). Absolute magnitudes of 3-vectors are denoted as k = |k|, and the 
unit vector in the direction of k is k = k/fc. 



II. DERIVATION OF THE PROPAGATOR FOR TRANSVERSE AND LONGITUDINAL GLUONS 



In this section, we derive the gluon propagator taking into account the fluctuations of the diquark condensate. A 
short version of this derivation can be found in Appendix C of Ref. || [see also the original Ref. j9)]. Nevertheless, 
for the sake of clarity and in order to make our presentation self-contained, we decide to present this once more in 
greater detail and in the notation of Ref. As this part is rather technical, the reader less interested in the details 
of the derivation should skip directly to our main result, Eqs. (|5^), (|57|), and (fsF 

We start with the grand partition function of QCD, 



Z = VA 



Z q [A] 



(2a) 



where 



Z q [A] 



T>ip Dip exp 



^7^ + /i 7 o+37^WV> 



(2b) 



is the grand partition function for massless quarks in the presence of a gluon field A^. In Eq. (|^), the space-time 

integration is defined as f = dr J v d 3 x , where V is the volume of the system, 7^ are the Dirac matrices, and 
T a = A a /2 are the generators of SU(N C ). For QCD, N c — 3, and A a are the Gell-Mann matrices. The quark fields 
ip are 4Af c Af/-component spinors, i.e., they carry Dirac indices a = 1, . . . , 4, fundamental color indices i = 1, . . . , N c , 
and flavor indices / = 1, . . . , Nf. The action for the gauge fields consists of three parts, 



Sa = S { 



S e [ + S 



FPG 



(3) 



where 



Sp2 



(4) 



is the gauge field part; here, F^ v = d^A^ — d u A a ^ + gf abc A b ^A c v is the field strength tensor, 
to gauge fixing, S s f , and to Fadeev-Popov ghosts, Sfpg: wu 1 be discussed later. 



The part corresponding 
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For fermions at finite chemical potential it is advantageous to introduce the charge-conjugate degrees of freedom 
explicitly. This restores the symmetry of the theory under fx — > — fx. Therefore, in Ref. 0], a kind of replica method 
was applied, in which one first artificially increases the number of quark species, and then replaces half of these 
species of quark fields by charge-conjugate quark fields. More precisely, first replace the quark partition function 
Z g [A] by Zm[A] = {Z q [A]} M , M being some large integer number. (Sending M — > 1 at the end of the calculation 
reproduces the original partition function.) Then, take M to be an even integer number, and replace the quark fields 
by charge-conjugate quark fields in M/2 of the factors Z q [A] in Zm[A]. This results in 



Z M [A] 



M/2 



exp 




<Z r (x)g- 1 (x,y)y r (y)+ / g9 r {x)A a Jx)t!£* r (x 



■r.y 



(5) 



Here, r labels the quark species and ^ r are 8iV c iV/-component Nambu-Gor'kov spinors 



WCr 



(lp r , IpCr) , 



(6) 



where ipcr = Ctff is the charge conjugate spinor and C = ij 2 jo is the charge conjugation matrix. The inverse of the 
8N c Nf x 8iV c iV/-dimensional Nambu-Gor'kov propagator for non-interacting quarks is defined as 



So" 1 



[Go"]" 1 
[Go]- 1 



where 



[G*]" 1 (x, y) ee -i (i^d* ± M7o) <$ (4) (x - y) 



(7) 



(8) 



is the inverse propagator for non-interacting quarks (upper sign) or charge conjugate quarks (lower sign), respectively. 
The Nambu-Gor'kov matrix vertex describing the interaction between quarks and gauge fields is defined as follows: 



where T£ = 7 T a and T£ = C{-)») T C- l T£ 



y"T„ T . 





o fa 



(9) 



Following Ref. |Q we now add the term J x ipcr{x) A + (x, y) ip r (y) and the corresponding charge-conjugate term 



f ip r {x) A (x,y)ipcr(y), where A = 7o(A + ) t 7o, to the argument of the exponent in Eq. (ph. This defines the 



quark (replica) partition function in the presence of the gluon field A% and the diquark source fields A + , A~ 

, M/2 ( M/2 f 



Z M [A, A+, A-] = / J] ™ r exp J ^ / * r (x) g-\x, y) V r (y) + f 

J r=1 r=1 Ux,y Jx 



g* r (x) A a Jx)T£ * r (x) 



(10) 



where 



[Go^]- 1 
A+ 



A- 
[Go-]- 1 



(11) 



is the inverse quasiparticle propagator. 

Inserting the partition function (|l(i| ) into Eq. (pa|), the (replica) QCD partition function is then computed in the 
presence of the (external) diquark source terms A^(x,y), Z — * Z[A + , A - ]. In principle, this is not the physically 
relevant quantity, from which one derives thermodynamic properties of the color superconductor. The diquark con- 
densate is not an external field, but assumes a nonzero value because of an intrinsic property of the system, namely 
the attractive gluon interaction in the color-antitriplet channel, which destabilizes the Fermi surface. 

The proper functional from which one derives thermodynamic functions is obtained by a Legendre transformation of 
lnZ[A + , A - ], in which the functional dependence on the diquark source term is replaced by that on the corresponding 
canonically conjugate variable, the diquark condensate. The Legendre-transformed functional is the effective action 
for the diquark condensate. If the latter is constant, the effective action is, up to a factor of V/T, identical to 
the effective potential. The effective potential is simply a function of the diquark condensate. Its explicit form for 
large-density QCD was derived in Ref. |Kj]. The value of this function at its maximum determines the pressure. The 
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maximum is determined by a Dyson-Schwinger equation for the diquark condensate, which is identical to the standard 
gap equation for the color-superconducting gap. It has been solved in the mean-field approximation in Refs. [fjl]-[l3"l| . 
In the mean-field approximation pi], 

A+(x, y ) ~ (ipcr(x)A(y)) , A~0,y) ~ (A(x)i>cr(y)) ■ (12) 

In this work, we are interested in the gluon propagator, and the derivation of the pressure via a Legendre transformation 
of lnZ[A + , A - ] is of no concern to us. In the following, we shall therefore continue to consider the partition function 
in the presence of (external) diquark source terms 

The diquark source terms in the quark (replica) partition function ( |Io|) could in principle be chosen differently for 
each quark species. This could be made explicit by giving A^ 1 a subscript r, A^ 1 — > A^. However, as we take the 
limit M — ► 1 at the end, it is not necessary to do so, as only A^ = A^ will survive anyway. In other words, we use 
the same diquark sources for all quark species. 

The next step is to explicitly investigate the fluctuations of the diquark condensate around its expectation value. 
These fluctuations correspond physically to the Nambu-Goldstone excitations (loosely termed "mesons" in the fol- 
lowing) in a color superconductor. As mentioned in the introduction, there are five such mesons in a two-flavor color 
superconductor, corresponding to the generators of SU(3) C which are broken in the color-superconducting phase. If 
the condensate is chosen to point in the (anti-) blue direction in fundamental color space, the broken generators are 
T4, . . . , T7 of the original SU (3) c group and the particular combination B = (1 + \/3T$)/ 3 of generators of the global 
U(1) B and local SU(3) C symmetry @. 

The effective action for the diquark condensate and, consequently, for the meson fields as fluctuations of the 
diquark condensate, is derived via a Legendre transformation of lnZ[A + , A - ]. In this work, we are concerned with 
the properties of the gluons and thus refrain from computing this effective action explicitly. Consequently, instead 
of considering the physical meson fields, we consider the variables in Z[A + , A - ], which correspond to these fields. 
These are the fluctuations of the diquark source terms A ± . We choose these fluctuations to be complex phase factors 
multiplying the magnitude of the source terms, 

A+(x,y) = V*{x)$+(x,y)V*(y) , (13a) 
A-{x,y) =V(x)$-(x,y)V T (y) , (13b) 



where 
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V{x) = exp 



<Pa(x)T a + ~J= (ps(x)B 



(14) 



The extra factor 1 / -\/3 in front of (p$ as compared to the treatment in Ref. |^] is chosen to simplify the notation in 
the following. 

Although the fields (p a are not the meson fields themselves, but external fields which, after a Legendre transformation 
of lnZ[A + , A - ], are replaced by the meson fields, we nevertheless (and somewhat imprecisely) refer to them as meson 
fields in the following. After having explicitly introduced the fluctuations of the diquark source terms in terms of 
phase factors, the functions < i> ± are only allowed to fluctuate in magnitude. For the sake of completeness, let us 
mention that one could again have introduced different fields tp ar for each replica r, but this is not really necessary, 
as we shall take the limit M — > 1 at the end of the calculation anyway. 

It is advantageous to also subject the quark fields ip r to a nonlinear transformation, introducing new fields \ r via 

Vv = VXr , Vv = Xr V f . (15) 

Since the meson fields are real- valued and the generators T4, . . . , T7 and B are hermitian, the (matrix- valued) operator 
V is unitary, V -1 = VL Therefore, the measure of the Grassmann integration over quark fields in Eq. ( |l0| ) remains 
unchanged. From Eq. ([l5]), the charge-conjugate fields transform as 

4>Cr = V* \Cr , -0CV = XCr V T , (16) 

The advantage of transforming the quark fields is that this preserves the simple structure of the terms coupling the 
quark fields to the diquark sources, 

^Cr(x) A + (x,y)tp r (y) = xcr(x)<S> + (x,y)xr(y) , tpr(x) A~{x,y)ipcr{y) = Xr(x)$~(x,y)xcr(y) ■ (17) 
In mean-field approximation, the diquark source terms are proportional to 
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® + (x,y) ~ (xcr(x)xr(y)) , $ (a:,y)~(xr(x)xcr(3/)) ■ (18) 

The transformation ( |l5|) has the following effect on the kinetic terms of the quarks and the term coupling quarks 
to gluons: 



Vv (i 7^ + (i 70 + g 7^ AJfT ) Vr = Xr (i y^/i + A* 7o + 7m w") Xr , (19a) 

V»cv (i 7^ - n 7o - 7m ^a T J) = Xcv (i 7^ - A* 7o + 7m w c) Xcv , (19b) 

where 

w" = V f (t + g A£T a ) V (20a) 

is the N c Nf x iV c 7V/-dimensional Maurer-Cartan one-form introduced in Ref. |l5[ and 

o£eeV t (id»-gA£T?) V* (20b) 

is its charge-conjugate version. Note that the partial derivative acts only on the phase factors V and V* on the right. 
Introducing the Nambu-Gor'kov spinors 

X r = ( ) , X r EE (Xr , XCr) (21) 



XCr 

and the 2N c Nf x 2iV c iV/-dimensional Maurer-Cartan one-form 



n"(a;,tf)=-i 

the quark (replica) partition function becomes 



Af/2 r M/2 

2i/[n, $+,$-] = / Y[vx r vx r explY, / 

J r=l [ r=l "^^'^ 



+ 7^(z, 2 /)] X r (y) } , (23) 



where 



* ' ( n- : ) • W 



o - 

We are interested in the properties of the gluons, and thus may integrate out the fermion fields. This integration 
can be performed analytically, with the result 

Z*[n, *+*-]= [det(5- 1 + 7A1 ^)] M/2 . (25) 

The determinant is to be taken over Nambu-Gor'kov, color, flavor, spin, and space-time indices. Finally, letting 
M — ► 1, we obtain the QCD partition function (in the presence of meson, ip a , and diquark, <J> ± , source fields) 



Z[<p,$ + ,$~] = J VA exp 



SA + iTrln(«S- 1 +7 At fi" 



(26) 

Remembering that is linear in A%, cf. Eq. (|22j) with (|20|), in order to derive the gluon propagator it is sufficient to 
expand the logarithm to second order in O^ 1 , 

i Tr In (S- 1 + 7|t fi") ~ 1 Tr ln^ 1 + ± Tr (5 7/1 fi") - ^Tr (5 7M fF 5 7 „fi'') 

= 5 [* + ,$-] + 5i[n,$ + ,$-] + 5 2 [n, *+,*-] , (27) 

with obvious definitions for the Si . The quasiparticle propagator is 

S=($.f-), (28) 



5 



with 

G ± = {[G^]- 1 -S ± } _1 , E ± = $^GJ$ ± , S ± = -GJ$ ± G ± . (29) 
To make further progress, we now expand and ojq to linear order in the meson fields, 

7 i 

W " ~ g A£ T a - ]T (0*V B ) T Q - -= (9^8) B , (30a) 

a=4 * 

~ -g A% Tj + J2 V**) T « + 7f BT ■ ( 30b ) 

a— 4 v 

The term Si in Eq. (^7| ) is simply a tadpole source term for the gluon fields. This term does not affect the gluon 
propagator, and thus can be ignored in the following. 

The quadratic term 5*2 represents the contribution of a fermion loop to the gluon self-energy. Its computation 
proceeds by first taking the trace over Nambu-Gor'kov space, 

$2 = -- / Tr Ci/)S [G + (x,y)-f fl uj fl (y)G + (y,x)j l/ uj''(x) + G~{x,y)j li u^(y)G~(y,x)'y v Wc } (x) 



+ E + (x,y) 7^(^)3 (y,x)j u uc(x) +5 (x, y) 7^(y) E + (y, x) i v uj u (x)\ . (31) 



The remaining trace runs only over color, flavor, and spin indices. Using translational invariance, the propagators 
and fields are now Fourier-transformed as 

G ± (x, y) = ^J2 e- iK<x ~ v) G ± (K) , (32a) 

K 

^(x,y) = ^^e- iK < x -^E ± (K) , (32b) 

K 

wHx) = J2 e ~ iPx wM ( p ) > ( 32 °) 
p 

u£(x)=J2e- iP - X "£(P) ■ (32d) 



Inserting this into Eq. (31), we arrive at Eq. (C16) of Ref. Q, which in our notation reads 

S 2 = -jE IV ',^ [ G+ ( K ) l^( p ) G + (K - P) lu w v {-P) + G-{K) 7 ^(P) G~(K - P) 7l *£(-P) 

+ ~ + (K) 7|tW "(P) S- (K - P) 7l/ w£(-P) + S- (K) 7„w&(P) S+(X - P) 7^"(-P)] . (33) 



4 

K,P 



The remainder of the calculation is straightforward, but somewhat tedious. First, insert the (Fourier-transform of 
the) linearized version ( |30| ) for the fields lo^ and This produces a plethora of terms which are second order in 
the gluon and meson fields, with coefficients that are traces over color, flavor, and spin. Next, perform the color and 
flavor traces in these coefficients. It turns out that some of them are identically zero, preventing the occurrence of 
terms which mix gluons of adjoint colors 1, 2, and 3 (the unbroken SU(2) C subgroup) among themselves and with 
the other gluon and meson fields. Furthermore, there are no terms mixing the meson fields ip a , a = 4, ... 7, with ipg. 
There are mixed terms between gluons and mesons with adjoint color indices 4, . . . , 7, and between the gluon field A^ 
and the meson field ipg. 

Some of the mixed terms (those which mix gluons and mesons of adjoint colors 4 and 5, as well as 6 and 7) can be 
eliminated via a unitary transformation analogous to the one employed in Ref. 0, Eq. (80). Introducing the tensors 

n£(p) = n^(p) = n^(p) = £■ | ]T Tr s [y* g+(k) y g+(k -p) + r g~(k) r g~{k p) 

K 

+ ^~-(K)Y^ + (K - P)+^~ + (K) 1 »~-(K - P)] , (34a) 

cf. Eq. (78a) of Ref. g, 



G 



II£ (P) = (P) = y ^ E Tr « [V G + W 7" G+ (K - P) + r G- (K) 7 " Go (K P)] 



2 V 



K 



cf. Eq. (83a) of Ref. §, 

n^(P) = II^(P) = ^E Tr s [V G+(K) ^ G+(K - P) + Y Go(K) ^ G~{K P)] 



2 V 



cf. Eq. (83b) of Ref. §, as well as 



K 



ngr(P) = |i^"(P) + |n^(P) , 

n 2 T 

IP" (p) = y_ _ ^ Tr s [y G+ (X") 7^ G+ (K - P) + 7^ G~ (if) 7" G~ (K - P) 
2 V K 

-^E-(K)^E + (K-P)- r S+ (K) Y E~(K- P)] , 
cf. Eq. (78c) of Ref. ](|, where ilp 1 " is the gluon self-energy in a dense, but normal-conducting system, 

K V iP) = yvE Tr * b" G » ^ l v G+{K -P) + r G - (K) Y Go (A" - P)] , 



(34b) 



(34c) 



(34d) 



(34e) 



(34f) 



K 



cf. Eq. (27b) of Ref. g], the final result can be written in the compact form (cf. Eq. (C19) of Ref. ||) 
1 V 



2 T 



EE 

P a=l 



A«(-P)--P /i ^(-P) 



IC(P) 



Al{P) + -P vV a {P) 
9 



(35) 



In deriving Eq. (|35J), we have made use of the transversality of the polarization tensor in the normal-conducting 
phase, Hq U (P)P v = P^Hq^(P) — 0. Note that the tensors for a = 1, 2, and 3 are also transverse, but those for 
a = 4, . . . , 8 are not. This can be seen explicitly from the expressions given in Ref. Q. The compact notation of Eq. 
(35) is made possible by the fact that ip a = for a — 1,2,3, and because we introduced the extra factor l/v3 in Eq. 
(14) as compared to Ref. ||. 

To make further progress, it is advantageous to tensor-decompose Various ways to do this are possible ||; 

here we follow the notation of Ref. B. First, define a projector onto the subspace parallel to P M , 



pp. pv 
—p2- 



E A " y 

Then choose a vector orthogonal to P M , for instance 

A^(^,^) = ( 5 --E-)/,, 
with / M = (0,p). Note that A^ 2 = —p 2 l p 2 /P 2 . Now define the projectors 

T>pv _ JT nn,v _ jy/j pv 1 pfi jy-v kuv _ uv _ -r>fJ.v _ 

N 2 ' ' 



(36) 



(37) 



(38) 



Using the explicit form of 7V M , one convinces oneself that the tensor A^ v projects onto the spatially transverse subspace 
orthogonal to P^, 



A 00 = A 0,: = , A ij = - (S ij - fpP) 



(39) 



(Reference || also uses the notation P£ v for A^ v .) Consequently, the tensor projects onto the spatially longitu- 
dinal subspace orthogonal to P' 1 , 



R oo P n oi PoP 



(40) 
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(Reference also employs the notation P£ v for B M ".) With these tensors, the gluon self-energy can be written in 
the form 



Km = Ka(P) + u h aa (p) b"" + ir aa (P) (T + n c aa (P) e 



(41) 



The polarization functions n^ a , n„ a , n„ „ , and II° Q can be computed by projecting the tensor onto the respective 
subspaces of the projectors (pq) and (J3q). Introducing the abbreviations 



Ka(P) = 2 



n^(p) , nl,(p) = ftii«(p)^ . 



■n^(p)A^ = -n* a (P), 



these functions read 

nS«(-P) 

Ka(P) 



K°a(P) 



2|n«(P)A 



2 TV 2 P 2 ^ Caw ~~ P 2 
1 



iC(p) 



pI+p 2 

PqP 



?lu e (P) 



n^(p)K + nL(P) 



rc (P) = ^ [ P g n°° (P) + 2 po p n« (p) ft + P 2 n£ a (p) 



(42) 

(43a) 
(43b) 

(43c) 
(43d) 



For the explicitly transverse tensor , the functions = nf x = 0. The same holds for the HDL polarization tensor 
n^. For the other gluon colors a = 4, . .. ,8, the functions n^ Q and n° a do not vanish. Note that the dimensions of 



and H% a are [MeV 2 ], while H^ a is dimensionless. 
Now let us define the functions 



n a n 



Al JP)=K/A a v {P) 



AMP) = 



P 2 



A a N {P) = 



Nt*A°(P) 
iV 2 



Note that Aff(-P) = -P'M£(-P)/P 2 , and A%(-P) = -N* A*(-P)/N 2 , since N» is odd under P 
fields ^4jj(P) and A^(P) are dimensionless. With the tensor decomposition ( fli] ) and the functions (J; 
becomes 



(44) 

-> -P. The 
1), Eq. @ 



S * = E E { ^1 n aa( P ) A ^ - A N(-P) U aa(P) N* A%(P) 

P a=l 



A\{-P) + -<f(-P) 

AU-P) + -ip a (-P) 
. 9 



n c aa (p) tv 2 p 2 a%(p) - a%(-p) n c aa (p) iv 2 p 2 



A? l (P) + -<p a (P) 



Ka(P)P 2 



Af(P) + -^(P) 



.9 



(45) 



In any spontaneously broken gauge theory, the excitations of the condensate mix with the gauge fields corresponding 
to the broken generators of the underlying gauge group. The mixing occurs in the components orthogonal to the 
spatially transverse degrees of freedom, i.e., for the spatially longitudinal fields, A a N , and the fields parallel to P M , 
^4jj . For the two-flavor color superconductor, these components mix with the meson fields for gluon colors 4, . . . , 8. 
The mixing is particularly evident in Eq. ( f45| ) . 

The terms mixing mesons and gauge fields can be eliminated by a suitable choice of gauge. The gauge to accomplish 
this goal is the 't Hooft gauge. The "unmixing" procedure of mesons and gauge fields consists of two steps. First, we 
eliminate the terms in Eq. ( f45[ ) which mix A^ and A?,. This is achieved by substituting 



IT (P)N 2 
A^P) = Ajf(P) + ™: > p) A%(P) . 



(46) 



(We do not perform this substitution for a = 1,2,3; for these gluon colors n^ Q = 0, such that there are no terms in 
Eq. ( |45[ ) which mix Ajj and A^)- This shift of the gauge field component A?, is completely innocuous for the following 

reasons. First, the Jacobian d(An , A^)/d(A« , A^) is unity, so the measure of the functional integral over gauge fields 
is not affected. Second, the only other term in the gauge field action, which is quadratic in the gauge fields and thus 
relevant for the derivation of the gluon propagator, is the free field action 
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'pi — 



1 V 



??EE^(- P ) (PV"-P^) A^(F)^-i^^^(-F)P 2 (A^+BH^(P) , (47) 



2 T 



P a=l 



2 T 



P a=l 



and it does not contain the parallel components A?< (P) . It is therefore also not affected by the shift of variables (|4(]) . 
After renaming Ajf — ► Ap the final result for 5*2 reads: 



1 V 

2 T 



E E { A ±. * (~ p ) n ~ ( p ) A i „ ( p ) - ^ ft L (^) ^ 2 ^ (^) 



P a=l 



. " 9 



Ka(P)P 2 



AUP) + ~<p a (P) 
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where we introduced 



n* a (p) = n£ a (p) 



K a (P)] 2 iV 2 P 2 



The 't Hooft gauge fixing term is now chosen to eliminate the mixing between A?i and ip a : 
1 V 



2A T 



EE 

P a=l 



p z Ai(-p)-\-n c aa {p)^{-p) 
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P*AUP)-\-W aa {P)v a {P) 



(48) 



(49) 



(50) 



This gauge condition is non-local in coordinate space, which seems peculiar, but poses no problem in momentum space. 
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Note that P 2 A?, (P) = P M A°(P). Therefore, in various limits the choice of gauge (|50| ) corresponds to covariant gauge, 



2A T 



(51) 



P a=l 



The first limit we consider is T, /i — + 0, i.e. the vacuum. Then, n® a = 0, and Eq. ( p0[ ) becomes (51). The second 
case is the limit of large 4-momenta, P — > oo. As shown in Ref. [Q, in this region of phase space the effects 
from a color-superconducting condensate on the gluon polarization tensor are negligible. In other words, the gluon 
polarization tensor approaches the HDL limit. The physical reason is that gluons with large momenta do not see 
quark Cooper pairs as composite objects, but resolve the individual color charges inside the pair. Consequently, 
I I" (P) P 2 — > P M IIq I/ (P) P v = for P — > oo and, for large P, the individual terms in the sum over P in Eqs. (|50| ) and 
( pl| ) agree. Finally, for gluon colors a = 1, 2, 3, Ii c aa = 0, since the self-energy H^" is transverse. Thus, for a = 1, 2, 3 
the terms in Eqs. ( |50t ) and (f5l| ) are identical. 

The decoupling of mesons and gluon degrees of freedom becomes obvious once we add ( |50| ) to (^) and (ft7|), 

4°2 + 5 2 + % = EE { [ p2 + Ka(p)] a"" ai.cp) 

P o=l 

-a^(-p) [p 2 + nL(P)l iv 2 A^(P) 
l 



Ajf(-P) 



A 
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^ + n c QQ (P) 



P 2 AMP) 



P 2 +Ka(P) 



Ka(P)<P a (P) 



Consequently, the inverse gluon propagator is 

a-^I(p) = [p 2 + Ka(P)} a^ + \p 2 + ii h aa (P) 



w 



lp2 + n c aa (p) 



Inverting this as discussed in Ref. one obtains the gluon propagator for gluons of color a, 

1 1 „„„ A 



K U a{P) = 



P 2 + U- a (P) 



A MV 



P 2 + ^ aa {P) 



P 2 + AIE aa (P) 



(52) 



(53) 



(54) 







For any A ^ 0, the gluon propagator contains unphysical contributions parallel to P M , which have to be cancelled by 
the corresponding Faddeev-Popov ghosts when computing physical observables. Only for A = these contributions 
vanish and the gluon propagator is explicitly transverse, i.e., P M A^(P) = A^(P)P„ = 0. Also, in this case the 
ghost propagator is independent of the chemical potential fi. The contribution of Fadeev-Popov ghosts to the gluon 
polarization tensor is then ~ g 2 T 2 and thus negligible at T = 0. We shall therefore focus on this particular choice 
for the gauge parameter in the following. Note that for A = 0, the inverse meson field propagator is 

D- 1 (P) = IF QQ (P) P 2 = P M UZ (P)P V , (55) 

and the dispersion relation for the mesons follows from the condition D~^(P) — 0, as demonstrated in Ref. Jl6| for a 
three-flavor color superconductor in the color-flavor-locked phase. 

The gluon propagator for transverse and longitudinal modes can now be read off Eq. (|4|) as coefficients of the 
corresponding tensors A MI/ (the projector onto the spatially transverse subspace orthogonal to P' 1 ) and B My (the 
projector onto the spatially longitudinal subspace orthogonal to P M ). For the transverse modes one has jl 

A «» (p) = P 2 + na a (P) = P 2 -n* a (P) ' (56) 



where we used Eq. (43a). Multiplying the coefficient of in Eq. (B4j) with the standard factor — P 2 /p 2 M, one 



obtains for the longitudinal modes 



P 2 1 1 



a°°(p) = -— ; = ; , (57) 

P 2 P 2 +Ka(P) P 2 -K°a(P) 



where the longitudinal gluon self-energy 



n° a °a(P)K a (P)-K a (P)Pi. 



TT UU ( P) — n 2 aa ' L aax L—l (£,0.) 

aa[ >- P p 2 W c P a {P) + 2p pX^ a {P)p l +p 2 Xii a {P) 1 ) 

follows from the definition of II^ , Eq. (|4S|), and the relations (|43]). The longitudinal gluon propagator must not 
be confused with the the 00-component of A^\ We deliberately use this (slightly ambiguous) notation to facilitate 
the comparison of our new and correct results with those of Ref. [Q, which were partially incorrect. The results of 
that paper were derived in Coulomb gauge, where the 00-component of the propagator is indeed identical to the 
longitudinal propagator (|57|). We were not able to find a 't Hooft gauge that converged to the Coulomb gauge in the 
various limits discussed above, and consequently had to base our discussion on the covariant gauge ( |5l| ) as limiting 
case of Eq. (|50|). 

To summarize this section, we have computed the gluon propagator for gluons in a two-flavor color superconductor. 
Due to condensation of quark Cooper pairs, the SU(3) C gauge symmetry is spontaneously broken to SU(2) C , leading 
to the appearance of five Nambu-Goldstone bosons. In general, these bosons mix with some components of the gauge 
fields corresponding to the broken generators. To "unmix" them we have used a form of 't Hooft gauge which smoothly 
converges to covariant gauge in the vacuum, as well as for large gluon momenta, and when the gluon polarization 
tensor is explicitly transverse. Finally, choosing the g aug e fixing parameter A = we derived the gluon propagator 
for transverse, Eq. (|5^), and longitudinal modes, Eq. (|5^ ) with (pq) . 



III. SPECTRAL PROPERTIES OF THE EIGHTH GLUON 



In this section, we explicitly compute the spectral properties of the eighth gluon. We shall confirm the results of 
Ref. 0| for the transverse mode and amend those for the longitudinal mode, which have not been correctly computed 
in Ref. 0. In particular, we shall show that the plasmon dispersion relation now has the correct behavior po — > m g 
as p — > 0. Furthermore, the longitudinal spectral density vanishes for gluon energies and momenta located on the 
dispersion branch of the Nambu-Goldstone bosons, i.e., for energies and momenta given by the roots of Eq. (^5|). For 
the eighth gluon, this condition can be written in the form P^ f[/" / (P) P v = ]9|Jl^|, since the HDL self-energy is 
transverse, P^ U^(P) P v = 0. 
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A. Polarization tensor 



We first compute the polarization tensor for the transverse and longitudinal components of the eighth gluon. To 
this end, it is convenient to rewrite the longitudinal gluon self-energy (pq) in the form 



n™(p)^nHp) + in 00 (P) 



n 00 (p) = p 2 



n 00 (p)W{P)- ir i (P)p, 



p^n 00 (p) + 2 PaP n 0i (p)p l +p 2 iL e (p) ' 



(59) 
(60) 



with n e {p) =pi w ] (p) pj . 

Let us now explicitly compute the polarization functions. As in Ref. we take T = 0, and we shall use the identity 



= V in S(x) , 



x + irj x 



(61) 



where V stands for the principal value description, in order to decompose the polarization tensor into real and 
imaginary parts. The imaginary parts can then be straightforwardly computed, while the real parts are computed 
from the dispersion integral 



1 [°° 
Ren(p ,p) = -V duj 

T J-oo 



Imil((jj, p) 
LO-po 



C 



(62a) 



where C is a (subtraction) constant. If ImII(a;,p) is an odd function of uj, ImII(— u>,p) = — Imn(cj,p), Eq. (|62af ) 
becomes Eq. (39) of Ref. @, 



^ /*oo 

ReU( PQ ,p) = -V 

7T 



dwIm±I odd (w,p) 



1 



s u)+Pq uj - Pq 

and if it is an even function of ui, ImII(— uj,p) = ImII(ix>,p), we have instead 

1 



lf°° ( 1 
Re±I(p ,p) = -V dui lmU even (u,p) 

T JO \U~ 



po u) + pa 



+ C , 



c 



(62b) 



(62c) 



Since the polarization tensor for the transverse gluon modes, II 



2 rr* _i_ I fr* 



88 — 3 n -r , 



II , has already been computed in 



Ref. 0, we just cite the results. The imaginary part of the transverse HDL polarization function reads (cf. Eq. (22b) 
of Ref. H) 



Imn*(P) = -vr i m- 



0{p-po) 



(63a) 



The corresponding real part is computed from Eq. (62b), with the result (cf. Eqs. (40b) and (41) of Ref. [Q) 



Kenl{P)= i -m 2 g 



Pa 



g+ll-g fin 
P V p J 2p 



Pq+P 



Po-P 



(63b) 



We have used the fact that the value of the subtraction constant is Cq = m 2 , which can be derived from comparing a 
direct calculation of Re H using Eq. (19b) of Ref. § with the above computation via the dispersion formula ( p2b| ). 
The imaginary part of the tensor II' is given by (cf . Eq. (36) of Ref. |7j ) 



Imn 4 (P) = -7T i m 2 q 9(pq -2^\ 9(E p - Po 



Pq 



1-^(1 + ^) Eft) - sM 1 - 2 ^ K(t) 



Pi 



0{po - E p ) 



Pi 



1-^(1 + s 2 )) E(a,t)-[1- 



Pq \ P 



P J Pq 



4 02 

P 2 Q-P 2 



-s z 1-2 



,Pl 



F(a,t) 



(64) 



where <fi is the value of the color-superconducting gap, E p = y/p 2 + 4</> 2 , t = yX — 40 2 /pq, s 2 = 1 — t 2 , a — 
arcsin[p/(tpo)], an d F(a,t), E(a,t) are elliptic integrals of the first and second kind, while K(t) = F(w/2,t) and 
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E(i) = E(ir/2,t) are the corresponding complete elliptic integrals. The real part is again computed from Eq. ( |62b] ). 
The integral has to be done numerically, see Appendix A of Ref. |?J for details. The subtraction constant is, for 
reasons discussed at length in Ref. [Q, identical to the one in the HDL limit, C* = Cq — m 2 . Finally, taking the linear 
combination n| 8 = § IIq + | IP completes the calculation of the transverse polarization function IT| 8 . 

In order to compute the polarization function for the longitudinal gluon, Ilgg , we have to know the functions IIq (P) , 
n 00 (P), n 0i (P)pi, and &(P). The first two functions, ng°(P) and fl 00 (P) have also been computed in Ref. 0. The 
imaginary part of the longitudinal HDL polarization function is (cf. Eq. (22a) of Ref. M ) 



ImU° °(P) = -nlm 2 g ^9(p-p 
z * p 



The real part is computed from Eq. (62b), with the result (cf. Eqs. (40a) and (41) of Ref. M) 



T 00/ 



Here, the subtraction constant is Cq° = 0. 



fl-^ln 


Po 


+ p 


) 


V 2p 


Po 


-p 





(65a) 



(65b) 



The imaginary part of the function n 00 is (cf. Eq. (35) of Ref. 0) 
Imn 00 (P) = -7T I m 2 9( Po - 2 4) ^ | 9(E P - Po ) E(i) + 0( Po - E p ) 



E(a,t) 



Po 



pI ~p 2 



(66) 



The real part is computed from Eq. ( 62b ), with the subtraction constant C 00 = Cg° = 0. Again, the integral has to 
be done numerically. 

It remains to compute the functions n 01 {P)pi and Il t (P). First, one performs the spin traces in Eq. ( |34c| ) to obtain 
Eqs. (102b) and (102c) of Ref. §. Then, taking T = 0, 



i 



i 



Po + ei + e 2 + irj po - ei - e 2 + if] 



H'GP) =~\ f J^js Yl [(l-eie 2 ki -k 2 ) +2 ei e 2 k! -pk 2 -p 



ei,e 2 =± 
X 



£l£2 - 66 



'J192 



2ei£ 2 



1 



(67a) 



(67b) 



where ki. 2 = k ± p/2, 



K Po + ei + £2 + i?l Po-ei-£2 + «^ / 
bfc. is the gap function for quasiparticles (ej = +1) or quasi- antiparticles (e^ = —1) with 



momentum k,, 6 = ejfcj — fi, and e l = yj (f + cj) 2 . 

One now repeats the steps discussed in detail in Section II. A of Ref. [0 to obtain (for p > 0) 



Im fl w (P) Vl = tt I m 2 g 9(p - 2 0) 4 { 0(E p - p ) [E(t) - s 2 K(t)] 



0(po - P P ) 



E(a,t) 



Po 



P 2 o-P 2 



s 2 F(a,t) 



(68a) 



ImII £ (P) = -n | m 2 0(p o - 2$ ^ { 0(£ p - Po ) [(1 + s 2 ) E(t) - 2 S 2 K(t)] 

+ 9(p - Ep) 



(l + s 2 )P(a,i)- — Jl- 



Po 



Po ~P 



-2s 2 F{a,t) 



One observes that in the limit — > 0, the functions ( pq ) approach the HDL result 

Imn?(P)K=^?m 2 4%-Po) , 



2 ' 9 p 2 

3 ™2P0 

P 3 



ImUi(P) = -7r-m 2 g ^9(p- Po ) 



(68b) 

(69a) 
(69b) 
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Applying Eq. ( pi] ) to Eqs. ( |67| ) we immediately see that the imaginary part of fl 0l (P)pi is even, while that of fr(P) 
is odd. Thus, in order to comp ute the real part of Tl 0l (P)pi, we have to use Eq. ( |62c| ), while the real part of IT(P) 
has to be computed from Eq. ( |62b| ) . When implementing the numerical procedure discussed in Appendix A of Ref. 
|jj for the integral in Eq. (62c), one has to modify Eq. (Al) of Ref. appropriately. 

Finally, one has to determine the values of the subtraction constants C° l and C e . We again use the fact that 
C°* = Cq 1 and C l = Cq, where the index "0" refers to the HDL limit. T he co rrespo ndin g constants are determined by 
first computing ReHQ l (P) fti and ReHp(P) from the dispersion formulas ( 62b ) a nd ( | 32c| ). The result of this calculation 
is then compared to that of a direct computation using, for instance, the result ( |65b| ) for Re I1q (P) and then inferring 
Ren^(P) pi and ReU e (P) from the transversality of lift". The result is C° l = = and C e = = m 2 g . 

At this point, we have determined all functions entering the transverse and longitudinal polarization functions for 
the eighth gluon. In Fig. [l] we show the imaginary parts and in Fig. |^ the real parts, for a fixed gluon momentum 
p = 4 0, as a function of gluon energy po (in units of 2 0). The units for the imaginary parts are — 3to^/2, and for 
the real parts +3m^/2. For comparison, in parts (a) and (g) of these figures, we show the results from Ref. for 
the longitudinal and transverse polarization function of the gluon with adjoint color 1. In parts (d), (e), and (f) the 
functions hi 00 , — Tl° l p*i, and IT are shown. According to Eq. (jfKj) these are required to determine II 00 , shown in part 
(b). Using Eq. (p9|), this result is then combined with the HDL polarization function Hq° to compute Ilgg, shown in 
part (c). Finally, the transverse polarization function for gluons of color 8 is shown in part (i). This function is given 
by the linear combination n 88 = | + -| n* of the transverse HDL polarization function and the function II' , 
both of which are shown in part (h). In all figures, the results for the two-flavor color superconductor are drawn as 
solid lines, while the dotted lines correspond to those in a normal conductor, 0^0 (the HDL limit). 

Note that parts (a), (d), (g), (h), and (i) of Figs. and agree with parts (a), (b), (d), (e), and (f) of Figs. 2 and 3 
of Ref. 0. The new results are parts (e) and (f) of Figs, [ljand ||, which are used to determine the functions in parts 
(b) and (c), the latter showing the correct longitudinal polarization function for the eighth gluon. In Ref. |t]], this 
function was not computed correctly, as the effect from the fluctuations of the condensate on the polarization tensor 
of the gluons was not taken into account. 




234012340123 
P /2* P /2* 

FIG. 1. Imaginary parts of the polarization tensors in a two-flavor color superconductor (solid lines) as a function of gluon 



energy p for fixed gluon momentum p = 4(f>. (a) Imllg, (b) Imfl 00 , (c) Imflgg, (d) Imfl 00 , 



(e) -Im^'ft, (f) Imtf, (g) 



Imlli!, (h) ImlT, (i) Imll| 8 . The corresponding results in the HDL limit, i.e. for = 0, are shown as dotted lines. 
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The singularity around a gluon energy somewhat smaller than po — 2 (f> visible in Figs. || (b) and (c) seems peculiar. 
It turns out that it arises due to a zero in the denominator of fl 00 in Eq. (|6(]), i.e., when P^ fl' 11 ' (P) P u — 0. As 
discussed above, this condition defines the dispersion branch of the Nambu-Goldstone excitations E™. Therefore, the 
singularity is tied to the existence of the Nambu-Goldstone excitations of the diquark condensate. 



p = 4(|) 




-2.0 1 ' 1 1 

0123401234012345 

P /24 P /2<» p /24 

FIG. 2. The same as in Fig. hi but for the real parts. 



B. Spectral densities 

Let us now determine the spectral densities for longitudinal and transverse modes, defined by (cf. Eq. (45) of Ref. 

0) 

Pss(Po,p) = -ImA^(po + "?,p) i P8s(Po,p) = -ImA| 8 (po + irj,p) (70) 

The longitudinal and transverse spectral densities for gluons of color 8 are shown in Figs. || (c) and (d), for fixed 
gluon momentum p = m g /2 and m g = 8 </>. For comparison, the corresponding spectral densities for gluons of color 1 
are shown in parts (a) and (b). Parts (a), (b), and (d) are identical to those of Fig. 6 of Ref. 0, part (c) is new and 
replaces Fig. 6 (c) of Ref. 0. One observes a peak in the spectral density around po — m g . This peak corresponds to 
the ordinary longitudinal gluon mode (the plasmon) present in a dense (or hot) medium. 

Note that the longitudinal spectral density for gluons of color 8 vanishes at an energy somewhat smaller than 
Po = m g/4- The reason is the singularity of the real part of the gluon self-energy seen in Figs. || (b) and (c). The 
location of this point is where P M n M "(P) P v = 0, i.e., on the dispersion branch of the Nambu-Goldstone excitations. 

Finally, we show in Fig. ^ the dispersion relations for all excitations, defined by the roots of 

p 2 -Ren°°Cp ,p) =0 (71a) 
for longitudinal gluons (cf. Eq. (47a) of Ref. Q|), and by the roots of 

pl-p 2 ~ReIll 8 (p ,p) =0 (71b) 
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for transverse gluons (cf. Eq. (47b) of Ref. 0). Let us mention that not all excitations found via Eqs. ( |7l| ) correspond 
to truly stable quasiparticles, i.e., the imaginary parts of the self-energies do not always vanish along the dispersion 
curves. Nevertheless, in that case Eqs. (|7l]) can still be used to identify peaks in the spectral densities, which 
correspond to unstable modes (which decay with a rate proportional to the width of the peak) . As long as the width 
of the peak (the decay rate of the quasiparticles) is small compared to its height, it makes sense to refer to these 
modes as quasiparticles. 
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FIG. 3. The longitudinal (a), (c) and transverse (b), (d) spectral densities for gluons of color 1 (a), (b) and 8 (c), (d). The 
gluon momentum is p — m g /2 and m g = 8<f>. For comparison, the dotted lines represent the corresponding HDL spectral 
densities. The poles of the spectral density corresponding to stable quasiparticles are made visible by using a numerically small 
but nonzero imaginary part. 
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FIG. 4. Dispersion relations for (a) longitudinal and (b) transverse modes for m g = 8 <j>. The solid lines are for gluons of color 
1, the dashed lines for gluons of color 8. The dotted lines correspond to the dispersion relations in the HDL limit. For both 
longitudinal and transverse gluons of color 8, the dispersion curves are indistinguishable from the HDL curves. The additional 
branch shown in (a) as dashed-dotted line is the one for the Nambu-Goldstone excitations, which appears as a zero in the 
longitudinal spectral density. 
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Fig. ^corresponds to Fig. 5 of Ref. 0. In fact, part (b) is identical in both figures. Fig. |J (a) differs from Fig. 5 (a) 
of Ref. [0, reflecting our new and correct results for the longitudinal gluon self-energy. In Fig. 5 (a) of Ref. 0, the 
dispersion curve for the longitudinal gluon of color 8 was seen to diverge for small gluon momenta. In Ref. it was 
argued that this behavior was due to neglecting the mesonic fluctuations of the diquark condensate. Indeed, properly 
accounting for these modes, we obtain a reasonable dispersion curve, approaching pq — m g as the momentum goes 
to zero. In Fig. || (a) we also show the dispersion branch for the Nambu-Goldstone excitations (dash-dotted). This 
is strictly speaking not given by a root of Eq. (fn]), but by the singularity of the real part of the longitudinal gluon 
self-energy. However, because this singularity involves a change of sign, a normal root-finding algorithm applied to 
Eq. ( |7l| ) will also locate this singularity. As expected [Q, the dispersion branch is linear, 

Po-^gP, (72) 
for small gluon momenta, and approaches the value po = 2 <fi for p — > oo. 



IV. CONCLUSIONS 



In cold, dense quark matter with Nf — 2 massless quark flavors, condensation of quark Cooper pairs spontaneously 
breaks the SU(3) C gauge symmetry to SU (2) c . This results in five Nambu-Goldstone excitations which mix with some 
of the components of the gluon fields corresponding to the broken generators. We have shown how to decouple them 
by a particular choice of 't Hooft gauge. The unphysical degrees of freedom in the gluon propagator can be eliminated 
by fixing the 't Hooft gauge parameter A = 0. In this way, we derived the propagator for transverse and longitudinal 
gluon modes in a two-flavor color superconductor accounting for the effect of the Nambu-Goldstone excitations. 

We then proceeded to explicitly compute the spectral properties of transverse and longitudinal gluons of adjoint 
color 8. The spectral density of the longitudinal mode now exhibits a well-behaved plasmon branch with the correct 
low-momentum limit p a — > m g . Moreover, the spectral density vanishes for gluon energies and momenta corresponding 
to the dispersion relation for Nambu-Goldstone excitations. We have thus amended and corrected previous results 
presented in Ref. jjj. 

Our results pose one final question: using the correct expression for the longitudinal self-energy of adjoint colors 
4, . . . , 8, do the values of the Debye masses derived in Ref. M change? The answer is "no" . In the limit po — 0, p — > 0, 
application of Eqs. (120), (124), and (129) of Ref. § to Eq. @ yields IT™ (0) = n°°(0), and the results of Ref. § 
for the Debye masses remain valid. 
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